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Abstract 

We conjecture that for a wide class of interacting particle systems evolving in dis- 
crete time, namely conservative cellular automata with piecewise linear flow diagram, 
relaxation to the limit set follows the same power law at critical points. We further 
describe the structure of the limit sets of such systems as unions of shifts of finite 
type. Relaxation to the equilibrium resembles ballistic annihilation, with "defects" 
propagating in opposite direction annihilating upon collision. 

1. Introduction 

Interacting particle systems evolving in discrete time have found many applications in recent 
years, especially in the field of road traffic modeling (see |I[ and references therein). While 
majority of traffic models based on interacting particle systems are stochastic, some features 
of the real traffic flow may be well described by purely deterministic models, as reported in 

B !!• 

An important characterization of a system of interacting particles is its fundamental 
diagram, i.e., the graph of the particle flux through a fixed point, in the equilibrium state, 
as a function of particle density. In many cases, such diagrams exhibit a discontinuity of the 
first derivative, and are very often piecewise linear. In the simplest cases there is only one 
such a discontinuity, and the fundamental diagram has the shape of an inverted "V" . Linear 
segments of the diagram can be interpreted as distinct "phases." 
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The discontinuity of the first derivative, to be further called a critical point, displays 
many features similar to critical phenomena known in statistical physics. The phenomenon 
of critical slowing down is especially apparent: when the density of particles approaches the 
critical density, the convergence to the equilibrium becomes slower. In what follows, we will 
present an evidence which strongly suggests that at the critical point the rate of convergence 
follows a power law with a universal exponent equal to — |. 



2. Interacting particle systems and cellular automata 

Discrete version of a totally asymmetric exclusion process is a simple yet often studied 
example of interacting particle system. |5|, 0, §, [J The dynamics of this process can be 
described as follows: particles reside on a one-dimensional lattice, with at most one particle 
per site. At each time t e N, each particle checks if the site to the right of its current 
position is empty, and if it is, it jumps to this site. The process is synchronous, meaning 
that all particles jump at the same time. 

Two alternative descriptions of this process are possible. We can label each particle with 
an integer n G Z, such that the closest particle to the right of particle n is labeled n + 1. If 
s(n, t) denotes the position of particle n at time t, the configuration of the particle system at 
time t is described by the increasing bisequence {s{n,t)} c ^ = _ OQ . The dynamics of the totally 
asymmetric exclusion process can be now stated as 

s(n, t + l) = s(n, t) + min-00 + 1, t) - s(n, t) -1,1}. (1) 

Another possible approach is to describe the process in the language of cellular automata. 
If lattice sites are labeled with consecutive integers ieZ, defining x(i,t) = 1 if the site i is 
occupied by a particle, and x(i, t) — if it is empty, the configuration of the particle system 
at time t is, in this case, described by the bisequence {x(i, t)}°Z_ OD . One can easily show 
that the dynamics of the aforementioned process is then given by 

x{i, t + 1) = / (x(i - 1, t),x(i, t),x{i + 1, t)) , (2) 

where / is defined by 

/(0, 0, 0) = 0, /(0, 0, 1) = 0, /(0, 1, 0) = 0, /(0, 1,1) = 1 (3) 
/(1,0,0) = 1, /(1,0,1) = 1, /(l, 1,0) = 0, /(l, 1, 1) = 1. 

The above definition can also be written in a more compact form as 

x(i, t + 1) = x(i, t) + min{x(z — 1, t), 1 — x(i, t)} — min{x(z, t), 1 — x(i + 1, t)}. (4) 

The evolution rule (£|) is often called cellular automaton rule 184, using the numbering 
scheme introduced by Wolfram [|TO|] . 

Consider now the general cellular automaton: 



x 



(i,t + 1) = f(x(i - r h t),x(i - n + l,t), . . . ,x(i + r r ,t)), (5) 
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where the function / : {0, l} r i+' v + 1 i— > {0, 1} is the evolution rule of the automaton, and the 
positive integers v\_ and r r are, respectively, the left and right radius of the rule. / will also 
be called a fc-input rule where k = r\ + r r + 1 is the number of arguments of /. A rule / 
is said to be conservative if for any periodic configuration of period L (i.e., a configuration 
such that x(i + L,t) = x(i) for every « G Z) we have 

L L 

£z(i,f + l) = ^z(i,t). (6) 

i=i i=i 

Every conservative cellular automaton rule can be viewed as the evolution rule of a 
system of interacting particles, just like rule 184 defined above, and a configuration of such a 
system can be represented by an increasing bisequence {s(n, t)}^!.^, where s(n,t) denotes 
the position of particle n at time t. A formal proof of this statement and an algorithm for 
construction of an analog of equation ([!]) for a given conservative rule / can be found in |[jfl . 
In [O] we have presented a survey of 4-input and 5-input conservative cellular automata 



rules and their properties. 

In order to describe the motion of the particles, we will use a visual representation of rule 
/ constructed as follows. List all relevant neighborhood configurations of a given particle 
represented by 1. Then, for each neighborhood, indicate the displacement of this particle 
drawing an arrow joining the initial and final positions of the particle. For rule 184, this 
would be 

10, 11, (7) 

where a circular arrow indicates that the particle does not move. The above notation is 
equivalent to saying that all particles which have empty sites immediately to their right 
jump one site to the right, while other particles do not move. 

We will normally list only neighborhoods resulting in particle motion, assuming that by 
default, in all other cases, particles do not move. Rule 184, for example, we will simply 

rv 

represented by 10. 

All particles do not necessarily move in the same direction. For example, there could 
be a process such that a particle will move one site to the right when the two nearest 
neighboring sites on its right are empty while, if its right neighboring site is occupied and 
its left neighboring site is empty, the particle moves one site to the left. In all other cases, 
the particle remains immobile. This rule can be visually represented by 

100, 011. (8) 

The main quantities of interest in this paper are the average particle velocity and flux at 
time t. For a periodic configuration of period L, the average velocity is defined as 



1 N 



n=l 



where N is the number of particles in a single period. The flux is defined as <f)(t) = pv(t), 
where p = N/L is the density of particles. In what follows, we assume that, at t = 0, 
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particles are randomly distributed on the lattice. When L — > oo, this corresponds to the 
Bernoulli product measure on Z, with lattice sites occupied by a particle with probability p, 
and empty with probability 1 — p. If b = &1&2 • • • is a block of finite length k, where, for 
all i e {1, 2, . . . , k} hi e {0, 1}, the probability of occurrence of block b in the configuration 
{x(i,t)}Z-oo wil1 be denoted by P t {b). 



3. Exact solution for deterministic traffic rule 



A natural extension of the totally asymmetric exclusion process with discrete time has been 
studied in connection with traffic models [I3|. Each site is either occupied by a particle, or 



empty. The velocity of each particle is an integer between and m. As before, if s(n,t) 
denotes the position of the nth car at time t, the position of the next car ahead at time t is 
s(n + 1, t). With this notation, the system evolves according to a synchronous rule given by 

s(n,t + 1) = s(n,t) +v(i,t), (10) 

where 

v(n, t) = min{s(n + 1, t) — s(n, t) — 1, m} (11) 

is the velocity of the particle n at time t. Since g = s(n+ 1, t) — s(n, t) — 1 is the gap (number 
of empty sites) between particles n and n + 1 at time t, one could say that each time step, 
each particle advances by g sites to the right if g < m, and by m sites if g > m. When 
m — 1, this model is equivalent to elementary cellular automaton rule 184 discussed in the 
previous section. 

If we start with a random initial (t = 0) configuration with particle density p, it is possible 
to obtain an exact expression for the flux at a later time t > 0, as shown in ||14|| : 

m = 1 - P - £ -i- ( (m + ^ 1 - p) m(m)+j - (12) 

j=1 -r \ + J 

In the limit t — > 00, we obtain the flux </>(oo) at equilibrium which is the piecewise linear 
function of p defined by 

0(oo) = (7 ifp<l/(m+l), 
y \ — p otherwise. 

We can see that for p < l/(m+ 1) the average velocity of particles at equilibrium is m, i.e., 
all particles are moving to the right with maximum speed m. The system is said to be in 
the a free-moving phase. When p > l/(m + 1), the speed of some particles is less than the 
maximum speed m. The system is in the so-called jammed phase. 

The transition from the free-moving phase to the jammed phase occurs at p = p c = 
l/(m + l) called the critical density. At p c , it is possible to obtain an asymptotic approxima- 
tion of (|T2|) by replacing the sum by an integral and using de Moivre-Laplace limit theorem, 
as done in |L4]]. At p c this procedure yields 



m ( (m + l) (m + l)t \ 

(00) - = J + ^ (e- ^ -e - ), (14) 
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Rule nr 


Motion representation 


3213933712 


rv rv rv 

0101 110 110 


3482385552 


rv rv rv 

0101 0110 110 


3486567632 


rv rx 

0101 110 


3163339916 


rv 

110 


3487613920 


rv rv 

0100 110 


3167653058 


•r> rv 

0001 1110 



Table 1: Examples of conservative rules with piecewise linear fundamental diagrams 



or, in other words, 0(oo) — <f){t) ~ t -1 / 2 . Here, by f(t) ~ git) we mean that lim^oo f(t)/g(t) 
exists and is different from 0. Power law convergence to the equilibrium at the critical point 
for rule 184 has been established in || and |I5" |. 

In the next section we will present some numerical results suggesting that this behavior 
is universal for a wide class of interacting particle systems evolving in discrete time. 



Numerical results for other CA rules 



Among conservative 5- input cellular automata rules investigated in ||12|| , a majority exhibit 
piecewise linear fundamental diagrams (i.e., graphs of 4>(oc) as a function of p). For the 
purpose of this article, we have selected six representative examples of such rules among the 
428 conservative 5-input rules. These six rules are given in Table 1, and their fundamental 
diagrams are shown in Figure 1. Points at which the slope of the fundamental diagram 
changes are referred to as critical points. In Figure 1, there are three rules with a single 
critical point, and whose fundamental diagrams are similar to rule 184. The remaining three 
rules have, respectively, 2,3, and 4 critical points. 

In order to illustrate that at a critical point, as a function of time t, the system converges 
to equilibrium as t~ 1 / 2 , we define the decay time as 



T = 

t=0 



$>(*)- (15) 



If the decay is of power-law type, the above sum diverges. 

For all rules in Table 1, we have performed computer simulations to estimate r. Results 
are shown in Figure 2, where r is plotted as a function of density for each rule. The value 
of r has been estimated by measuring <j)(t) for t — 0,1, ... , 1000, and truncating the sum 



(|15| ) at t = 1000. Comparing Figures 1 and 2 we clearly see that r diverges at critical 
points. In order to verify that indeed |0(oo) — <p(t)\ ~ t~ x l 2 at critical points, we have 
plotted (for all critical points) the numerical values of |0(oo) — <p{t)\ as a function of t using 
logarithmic coordinates. Figure 3 shows an example of such a graph, for the critical point of 
the first rule in Table 1 (code number: 3213933712). The (p, 0) coordinates of this critical 
point in the fundamental diagram are (|, |). Assuming that |0(oo) — 4>(t)\ ~ t~ a , we have 
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Figure 1: Fundamental diagrams for rules from Table 1. 
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Figure 2: Decay time as a function of density for rules from Table 1. 
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Figure 3: Logarithmic plot of |0(oo) — <f){t)\ as a function of t for rule 3213933712. Data 
points (+) represent computer simulations, while the dashed line represents the best fit 
(slope = -0.501 ±0.01) 

determined the exponent a as the slope of the straight line which best fit the logarithmic 
plot of \4>(oo) — <f)(t)\. Table 2 shows our results for all the critical points of all the six 
rules of Table 1. All these exponents are remarkably close to |, suggesting that t _1//2 might 
be a universal law governing the rate of convergence to equilibrium at critical points of 
conservative cellular automata with piecewise linear fundamental diagrams. 



5. Shifts of finite type 

In order to understand why the decay to equilibrium, or, in other words, approach to the 
limit set follows the same law for different rules, we now describe limit sets of conservative 
rules using the concept of a shift of finite type (SFT). In symbolic dynamics and coding 
theory ||16|| , SFT is defined as follows: let B be a finite set of finite blocks (consecutive sites), 
e.g. B = {00, 11, 101}. Consider a set of all bi-infinite configurations in which blocks of the 
set B do not appear ("forbidden blocks"). Such a set, denoted by SFT(S), is called a shift 
of finite type. 

By analyzing frequencies of occurrences of finite blocks in conservative CA, one can ob- 
serve that limit sets of many conservative CA are unions of two or more SFT's. Each of 
these SFT's correspond to a distinct "phase", or a straight line segment in the fundamen- 
tal diagram. Critical points correspond to a set of configurations common to two SFT's. 
Forbidden block sets of each phase completely determine the fundamental diagram. 

As an example, consider the deterministic traffic rule defined by the relations (|10|) and 
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Rule nr 


Critical point (p, </>) 


a 


3213933712 


(2/3,2/3) 


0.501 


3482385552 


(1/2, 1/2) 


0.517 


3486567632 


(2/3, 1/3) 


0.507 


3163339916 


(1/2,0) 


0.506 


(2/3, 1/3) 


0.504 


3487613920 


(1/3, 1/3) 


0.502 


(1/2,0) 


0.539 


(2/3, 1/3) 


0.504 


3167653058 


(1/4, -1/4) 


0.486 


(1/3,0) 


0.514 


(2/3,0) 


0.514 


(3/4, 1/4) 


0.486 



Table 2: Values of the exponent a at critical points of rules from Table 1. Simulations 
performed using a lattice size equal to 5 x 10 5 . 



(|TT|), assuming m — 2. We mentioned that in the free-moving phase, all particles move 
with maximum speed m — 2, which implies that l's are always separated by two or more 
0's. Hence, the blocks 11 and 101 cannot be found in a configuration belonging to the 
limit set of the free-moving phase. In the jammed phase, on the other hand, blocks of 
zeros of length 3 or more cannot occur. The limit sets of these two phases are therefore, 
respectively, SFT(11, 101) and SFT(000). Note that these two sets are not completely dis- 
joint, but have three common elements, namely, ...100100100100..., ...001001001001..., and 
...010010010010... 

Consider now the phase PI = SFT(11,101). It is clear that, for any configuration in 
PI, p < 1/3, since two ones must be always separated by two or more zeros. We know 
that the flow for the above rule is defined as <fi(t) = 2P(100) + P(101). We also know that 
P(ll) = P(101) = 0. The consistency condition gives P(100)+P(101)+P(110)+P(lll) = p. 
Hence, P(100) = p, and, in this phase, 0(oo) = 2p. 

In phase P2 = SFT(000), P(000) = 0, and one can show that 0(oo) = 1 - p - P(000). 
The proof can be found in [[14j] (it involves manipulations of the consistency conditions only). 
Hence, = 1 — p. In P2, the largest number of zeros separating two ones is 2, hence the 
minimum value of the density p is |, so 0(oo) = 1 — p is valid for all p > |. 

Note that we derived the fundamental diagram given by equation (|13|) from the struc- 
ture of the limit set only. This can be done for all conservative CA with piecewise linear 
fundamental diagrams. As mentioned earlier, different straight line segments of the funda- 
mental diagram correspond to different SFT components of the limit set, while configurations 
common to two SFT's define critical points. 

At the critical point, the limit set consists of a finite number of configurations, yet there 
is an infinite number of initial configurations with a given density. Time evolution of the 
system governed by a conservative rule can, therefore, be viewed as a transition form an 
infinite configuration space to a finite configuration space. While details of this transition 
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are different for different rules, it appears that the rate of convergence follows always the 
same law, gued in the previous section. 

In order to test this interpretation of the dynamics of conservative CA, we have deter- 
mined the spatial measure entropy of configurations at time t. The spatial measure entropy 
of a limit set configuration is defined as 

si * k) = -\ E p ~ (6) log Poo (6), (16) 
£>e{o,i} fe 

where the sum runs over all blocks of length k, and Poo(b) denotes the probability of oc- 
currence of block b in the limit set. If, at a given p, the number of configurations in the 
limit set is finite, s(k) should go to zero with k — > oo. Figure 4 shows s(10) plotted as a 
function of the density p for all the six rules of Table 1. Comparing Figures 1 and 4, one can 
easily notice that, at critical points of fundamental diagrams, s(10) takes minimum values, 
confirming our observations regarding the structure of limit sets. 



6. Annihilation of defects 

The dynamics of rule 184 can be described in terms of soliton-like localized structures propa- 
gating in periodic background |T7[ More formally, it can be shown that rule 184 is equivalent 



to ballistic annihilation process [pTSfl . This equivalence has been used in ||, and investigated 
in detail in JT5[]. The ballistic annihilation process involves two types of particles moving, on 
a line, at constant speed in opposite directions. When two particles collide, they annihilate. 
Using the central limit theorem, Elskens and Frisch fl8|] demonstrated that if the initial 
distribution of particles is balanced, then, the fraction S(t) of surviving particles at time t 
behaves as t~ x l 2 . 

For rule 184, the "background" on which the defects are moving are periodic configura- 
tions of alternating zeros and ones, . . . 01010101 .... If we call A-type defect consecutive l's, 
and -B-type defect two consecutive O's, it is easy to show |1 that, for each iteration of rule 
184, defects type A move to the left, while defects of type B move to the right. When they 
collide, the reaction A + B — ► background, i.e., annihilation, takes place. 

A very similar construction can be made for the process defined by equation (|10D and 



([Til) • Here, the background consists of ones separated by m zeros. We can also define two 



types of defects, but with an additional subscript identifying their length. If two consecutive 
l's are separated by a cluster of zeros of length m + k {k > 0), this is the defect of type 
B, denoted B^. Similarly, if two l's are separated by a cluster of zeros of length m — I 
(m > I > 0), such a cluster constitutes an A- type defect, denoted A t . A-type defects move 
to the left with speed 1 while -B-type defects move to the right with speed m. A collision 
between defects A\ and B^ results in a defect Ai_^ if I > k, and Bk_i if k > I. If k = I, the 
two defects annihilate. Figure 5 shows an example of such collisions. 

The dynamics of other conservative rules with piecewise linear fundamental diagrams can 
also be interpreted as interaction of defects propagating in a periodic background. However, 
the number of different types of "defects" and "backgrounds" is typically larger than in the 
simple cases described above. This can be seen in Figure 6, which shows spatiotemporal 
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Figure 4: Spatial measure entropy for block length k = 10 as a function of density for rules 
of Table I. 
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Figure 5: Collision of defects in deterministic traffic rule defined by relations ( [TOD and (|TT|) 
with m = 2. Black squares represent lattice sites occupied by a particle. First, A\ collides 
with B 3 , resulting in B 2 . Several iterations later, B 2 collides with A 2 , and both defects 
disappear. 

diagrams for rule 31633399164 (110) at p — ^ and p = |. "Defects" propagating in opposite 
directions and annihilating upon collision can be clearly identified in both diagrams. At p — \ 
all defects eventually disappear resulting in a periodic configuration . . . 01010101 .... Simi- 
larly, at p = | annihilation of all defects results in a periodic configuration . . . 110110110 .... 
Since these processes can be considered a generalization of a simple ballistic annihilation, 
one can conjecture that their rate of convergence to equilibrium should be the same as for 
the ballistic annihilation, in agreement with the results of computer simulations presented 
in a previous section. 

7. Conclusion 

We have presented numerical evidence suggesting that, for conservative cellular automata 
with piecewise linear flow diagrams, relaxation to equilibrium at a critical point follows a 
universal power law with exponent — |. The universality of this critical behavior is related 
to the structure of limit sets of such rules, which can be described as unions of shifts of finite 
type. At critical points the spatial measure entropy of the limit set goes to zero, which means 
that the limit set consists of a finite number of configurations. Moreover, the dynamics of 
such rules can be viewed as interaction of "defects" propagating in opposite directions and 
annihilating upon collision, just like in a simple ballistic annihilation process, for which the 
power law t^ 1 ^ 2 behavior has been established analytically. 
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